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Abstract 

Wc pose and solve a problem concerning consistent assignment of 
quantum probabilities to a set of bases associated with maximal pro- 
jective measurements. We show that our solution is optimal. We also 
consider some consequences of the main theorem in the paper in con- 
junction with Gleason's theorem. Some potential applications to state 
tomography and probabilistic quantum secret-sharing scheme are dis- 
cussed. 



1 Assignment of quantum probabilities 

Consider a finite quantum system 5f of dimension n. Define a quantum prob- 
ability assignment (QPA) for 5? to be a set & = {{Bi, Si), ... , (B m , S m )} 
consisting of m pairs (Bk, Sk) where 



1. each Bk = (Pj , ■ ■ ■ , P*) is a sequence of n 1-dimensional projection op- 

5^ I \ I / . / \-y-y / \: t\ i in rliTi/T +■ /-\ Lit, / \ v J- 1 1 ^ , v i / , v v v i ol K^cin I ^r^\ I /'V \ 



erators Pf = \af) (oq\ corresponding to an orthonormal basis \af) , . . . , |a* 
and 



2. each Sk = {pi,P2, ■ ■ ■ ,Pn) ^ s a sequence of n real numbers such that 
p\ > and Y.iP'i = !■ 

Clearly, the sequences Sk can be considered as distributions on a probability 
space with n elements (atomic events). 
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Suppose now it is claimed that these probabilities correspond to pro- 
jective measurements in these m bases. Thus, the claim is that there is a 
(mixed) quantum state p (a positive semidefinite matrix with trace 1) such 
that 

Tx(pP?) = for all 1 < k < m, 1 < i < n (1) 

If such a state exists we say that the collection is a consistent quantum 
probability assignment. How do we verify this claim? Potentially m could 
be infinite. A related question is the following 

Problem. What is the smallest number r such that if all the probability 
assignments ^ C & with < r i\ c S\ denotes cardinality of £f ) are consis- 
tent, then the whole collection & is consistent? 

The number depends on the assignment & . So we ask what is the largest 
value of r(j^~) as we vary J^~. This value depends only on the dimension n. 
Call it the consistency number r n . Equivalently, r n is the smallest number 
such that for every collection of quantum probability assignments J^, if every 
subset of & of size r n is consistent then & is consistent. This problem 
first came up in our attempt at formal axiomatization of a logic of finite- 
dimensional quantum systems pp. The problem we address in this paper is 
the calculation of the consistency number. 

We mention that a related problem was posed by Bell and Kochen and 
Specker [21 [3] in their analysis of non-contextual hidden variable theories. 
Projection operators of rank 1 are the quantum analogues of classical prepo- 
sitional variables. The latter can take two values, say and 1, corresponding 
to true and false. The question then is whether there is a consistent prob- 
ability assignment to the whole collection such that the p\ take value 1 or 
0. Clearly quantum mechanics prohibits this in case of sets of incompatible 
projection operators. The question posed by these authors is whether such 
assignment is possible for an arbitrary collection of bases in some hidden 
variable theory. 

The paper is structured as follows. In Section [2] and El we prove that 
r2 = 4 and r n = n 2 — n + 1 for n > 3. We first prove a weaker result that 
if every subset consisting of r n + 1 = n 2 — n + 2 probability assignments 
is consistent then so is the whole set & . The proof in the case of r n is 
surprisingly much harder. There is one exceptional case where one has to 
investigate in detail the structure of the bases themselves (unlike the case 
for r n + 1). But we get extra dividends in gaining information about the 
measurement bases. In fact, there is essentially one case where we need the 
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maximal value n 2 — n + 1. We use this in Section [3] to construct examples to 
demonstrate that the number r n is indeed optimal. More precisely, if r' < r n 
there exist quantum probability assignments & such that all subcollections 
of & of size r' are consistent but & itself is not consistent. The tricky part 
in the construction of these examples is that we have to ensure that positive 
definite matrices with trace 1 (states) exist which satisfy the assignments 
corresponding to every subcollection of size r' . We conclude the section with 
a simple application of the main result in conjunction with Gleason's famous 
theorem [BJ- In Section [3] we outline some applications of our constructions 
to (probabilistic) secret sharing and state tomography. Section [5] makes 
some concluding remarks. 

2 Solution to the problem 

This section is devoted to establishing upper bounds for the consistency 
number. We assume first that the number of assignments m is finite, but 
we will show later that this assumption can be dropped. More precisely, we 
prove: 

Theorem 1. Let H be the Hilbert space of an n- dimensional quantum sys- 
tem. Define a quantum probability assignment & = {(Si, Si), ... , (B m , S m )} 
consisting of m pairs {Bk,Sk) where 

i. Bk = (Pi, ■ ■ ■ , Pn) is a sequence ofn projection operators Pf = \cq) (a^ 
corresponding to an orthonormal basis \a\) , . . . , \cQ), and 

ii. Sfc = (Pi,P2, ■ ■ ■ ,Pn) i s a sequence of n real numbers such that p\ > 



Then & is a consistent quantum probability assignment to an n-dimensional 
quantum system if and only if every subset C & of size at most 



Note that arbitrary fe-subsets may be separately consistent. How- 
ever, the proof below will produce a state in which the whole of & is a 
consistent assignment provided the conditions of the theorem are satisfied. 
Observe also that the problem of assignment arises because of the unitary 
transformation connecting two bases. We will later see that these bounds 



and J2iPi = 1- 





is consistent. That is, r n < R 
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are optimal in the sense that for a positive integer r' < R n there are bases 
and corresponding probability assignments such that every subset of size r' 
is consistent but the whole set is inconsistent. 
We will start with a relatively simpler result. 

Lemma 1. With the notation as above, & = {(B\, S±), . . . , (B m , S m }} is 
a consistent probability assignment to an n- dimensional quantum system if 
and only if every subset W C & of size at most n 2 — n + 2 is consistent. 

Proof. Let & = {(B\, S\), . . . , (B m , S m )} where Bi are bases and Si the 
corresponding probability assignments. We may assume that m > n 2 — n + 2. 
Let I denote the unit matrix of order n. Now for each 1 < k < m we have 
a set of consistency equations 

Tr(pl) = Tr(p) = 1 
1r(pP t k ) = >p\ for 1 < i < n - 1 

in the real variables xjk and yjk where pjk = Xjk + iyjk, the (jk) th entry 
of the quantum state p. Note that we omit the equation Tr (pP„) = p k for 
the case i = n since this follows from others usine; the facts that ^ Pf = I 
and Y2iPi = 1- The hermiticity of p implies that there are n 2 independent 
parameters. Any one set of equations corresponding to a single basis assign- 
ment fixes the trace condition. We call the set of equations corresponding 
to the assignment to the kth. basis the feth block, P>k- So we have n 2 in- 
dependent real variables and any block, say the first, has n independent 
equations. That is, the rank of the coefficient matrix of the first block is n. 
But after the normalization (Tr(/j) = 1) is fixed, each block from the second 
block onwards can contribute at most n — 1 to the overall rank. Given an 
arbitrary set A of assignments let rank(^4) denote the rank of the coeffi- 
cient matrix of the corresponding set of linear equations ([3]) . Since p has n 2 
real variables rank(^l) < n 2 . Now consider a subset £f C & of cardinality 
/ > n 2 — n + 2. Call such subsets Z-subsets of & . We will show that if every 
/-subset is consistent then any (I + l)-subset is consistent. Then induction 
will complete the proof. Hence, we have to show that given any (I + 1)- 
subset = {(Bi, Si) : 1 < i < I + 1} it is consistent if every /-subset of 
is consistent. 

By assumption ^ = — {(Bi, Si)} is consistent for all 1 < i < I + 1. 
Now consider the sets Sfy = % — {(Bj, Sj)}. Suppose rank(5fy) < n 2 . Since 
1 5% I > n 2 — n + 1 this implies that at least one of the blocks say Bk in Sf^ is 
dependent on the rest in the sense that rank(?fjj — {{Bk, Sk)}) = rank(^). 
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Suppose this is not the case. Then we must have 

n = rank({(Bi, Si)}) < rank((#i, Si), (B 2 , S 2 ) })<••• < rank(%) < n 2 

But this is impossible since > n 2 — n + 1. We have shown that there 
is some k such that (B k ,S k ) £ S% and rank(^j — {(£>fc,Sfc)}) = rank(^j). 
Let = @ij — (B k ,S k ). Then every equation in the block B k can be 
written as a linear combination of equations of Hence consistency of 
&ij implies that any solution to ^ is a solution for ^ . The cardinality of 
&ljU{{(Bi, Si) , (Bi, Si)} is / and by hypothesis it is consistent. But from 
what we have shown above a solution p to this system is also a solution to 

U{{(^i> (-^j) (-^fc) Sfc)} = The proof is complete in this case. 

Next consider the second alternative, rank(^j) = n 2 . This is the maxi- 
mal rank since there are n 2 variables. Hence, any solution to the system £fy 
is unique. This implies that ^ and &j have the same unique solution. But 
then so does o€ and the proof is complete. □ 

The following lemma provides one of the key elements in improving the 
bound so as to obtain our main result. 

Lemma 2. Let B = {ei, . . . , e n } and B' = {(3i, . . . , /3 n } be two bases in an 
n- dimensional Hilbert space. Let B and B' be, respectively, the associated 
matrices of coefficients for the systems of real equations 

Trflej) (ei\p) = p h and Tr(|/3j) p) = q h (both for i = 1, ... ,n) (4) 

where Y2i Pi = Si Qi = 1 an d where p is positive semidefinite matrix and 
the entries of p are treated as unknown variables. Suppose the rank of the 
combined system B and B' is n + 1, exactly 1 higher than the rank of B 
(or B'). Then two of the vectors from B' , say, (3\ and $2 He in the plane 
determined by two vectors from B. Assuming (without loss of generality) that 
the later is spanned by e± and €2 we also have (ej| /3j) 7^ for i,j 6 {1,2}. 
Further, the set {f3j : 3 < j < n} is a permutation of {ej : 3 < j < n}. 

Proof. Writing the matrix elements of p in the B basis we observe that the 
equations corresponding to B determine the diagonal elements pa along with 
the constraint Yli Pa = 1- Now let = Yli <kj£i- Then the equations in B' 
are equivalent to the following. 

\cji\ 2 Xjj + ^2(cjic M + Cjic k i)x jk + i(cjic k i - Cjic k i)y jk = q t where 

3 3<k 

Pjk = (ej\p\ek) = Xjk + Wjk 
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Using B we can eliminate the diagonal terms from these equations (corre- 
sponding to j = k). Then the hypothesis that the rank of the combined 
system B and B' is exactly one more than that of B alone implies that 

CjlCkl = aimCjmCkm Vj / k (5) 

Here ai m are real constants depending on the basis vectors /?/ and /3 m or 
equivalently the Zth and mth row vectors of B' . Now, among the fas there 
must be at least two vectors which have at least two nonzero coefficients 
when expressed in the vectors of B. Otherwise, B and B' are the same basis 
apart from permutation. We may take these two vector to be fa and fa. 
We may also assume without loss of generality that en and C21 are both 
nonzero. Suppose for fa the coefficients Cj2 and 0^2 are nonzero. Then, 
addition of either of the equations corresponding to fa or fa to the system 
corresponding to B increases rank by 1. The equation [5] then shows that 
a\2 and a2i(= 01^2) are both nonzero. Hence, C21 and C22 must be nonzero. 
Then equation ([5]) implies that C12 and C22 are both nonzero. Using equation 
[5] again we infer that there is a nonzero constant 7 such that Cj2 = JCji for 
j = 3, ... ,n. If a third basis vector, say fa, had two nonzero coefficients 
a similar type of argument would show the existence of a nonzero constant 
7' such that Cj3 = 7'cji, j = 2,3, ... ,n. But then fa, fa and fa will not 
be linearly independent. Hence, all fa, i = 3, ...,n must have exactly 
one nonzero coefficient cti- Since they are unit vectors this must 1 (apart 
from an inconsequential factor of modulus 1). That is, fa = ej r Since 
the fa's constitute an orthogonal basis ej i cannot be any of the whose 
coefficients Cki (and 0^2) is nonzero. In particular, ji 7^ 1 or 2. But there 
are n — 2 such /3j's. Hence, the set {fa, fa, ■ ■ ■ , fa} must be a permutation of 
{€3, €4, . . . , e n } and {fa, fa} form an orthonormal basis in the 2-dimensional 
space spanned by ei and 62- Moreover, from the fact that c\\, C21, C12 and C22 
are all nonzero we conclude that the (ej| fa) 7^ for i,j G {1,2}. The proof 
of the lemma is complete. □ 

We now improve Lemma Q] so as to yield a proof of Theorem (U by 
considering several special cases. 

Proof. The implication from left to right is trivial — any subset of a consis- 
tent set is consistent. To prove the converse, we may assume m > R n . 

The case n = 2 is covered by Lemma [TJ So we assume that n > 3 
and that any i? n -subset of is consistent. The problem can be restated 
as follows. Find a positive semidefinite matrix p satisfying the set of linear 
equations ([3]). There are n 2 independent real variables corresponding to the 
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real and imaginary parts of entries of p and hence the real rank of system 
is < n 2 . Any one of the blocks corresponding to one basis assignment has 
rank n. For each of the rest of blocks of equations corresponding to the pair 
{Bk,Sk) we need to verify n — 1 equations in each block since the last one 
is guaranteed by the condition Yl,iPi = 1- Again using Lemma [T] it suffices 
to show that every (R n + l)-subset of & is consistent. Assume otherwise: 
there is a set = {(£>i, Si), . . . , (Br„+i, Sji n+ i)} such that £f is inconsistent 
while every proper subset of is consistent. 

Recall that the rank(^) of a quantum probability assignment is defined 
to be the rank of the coefficient matrix of the consistency equations fl3J) 
corresponding to <S . The set 5f' = {(B±, Si), . . . , (Br^, SrJ)} is consistent 
by hypothesis. Suppose that for some pair (Bi,Si), (Bj,Sj) G Sf', the rank 
of the combined system satisfies 

rank({B 4 , Si), (B^Sj)}) > n + 2 . (6) 

Consider the chain of inequalities 

rank({(^,S 4 )}) = n < n + 2 < rank((^,S 4 ), {Bj,Sj)}) < 
rank({(^, Si), (Bj, Sj), (Si, Si)}) < . . . < rank(^) < n 2 

where in each term we add exactly one new basis assignment to the previous 
set. Since there are R n = n 2 — n + 1 terms at least one of the relations < 
must be an equality. Hence, the equations corresponding to at least one of 
the assignments, say, (Bk,Sk) must be dependent on the rest. Thus, every 
solution of &'-{(Bk, S k )} must be a solution of Sf'. Now, $ k = Sf-{(B fc , S k )} 
is consistent being a set of cardinality R n . From what we have just seen 
any solution to % is also a solution for (Bk,Sk)- Hence §f is consistent, a 
contradiction. We observe that in general we arrive at the same conclusion 
if any of the relations in the chain 

ranktt^,^)}) < rank({(B il , S h ), (B i2 , S h )}) < ... < rankgf (7) 

where (ii,t2, • • • , «R n ) is a permutation of (1,2,..., R n ) is an equality. Note 
also that if rank(^') < n 2 then an argument similar to the above would 
show that some assignment {B^, S&} must be dependent on the rest without 
requiring the condition in (|6|). We can then show that & is consistent. 
Hence, to complete the proof we have to consider only the following case. 

Condition 1. For every distinct pair in 1 . . . R n + 1, 

rank({(^, Si), (Bj,Sj)}) = rank({(Bi, Si)}) + 1 = n + 1 (8) 
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and for every i? n -subset <S' = {(5^,5^), (B i2 ,S i2 ), . . . (B iRn ,S iRn )} of <S , in 
the chain 

rank({(£ n ,S n >}) < rank({(B il , S h ), (B i2 , S i2 )}) <■■■< rank(Sf) (9) 
the rank increases exactly by 1 at each step. 

This condition implies that the set of consistency equations correspond- 
ing to every itl n -subset Sf' of have rank n 2 . Let us write the basis vectors 
of the basis B^ as 

B k = {$,$,...,/3*} 

Without loss of generality we may take the basis B\ to be any fixed basis, say, 
the "computational" basis {ei, . . . ,e n }0 Then the relation (JSj) and lemma 
[2] imply that for the pair of distinct indices (1, j) from l..R n + 1 there exist 
pairs of distinct indices (r, s) and (j r ,js) from l..n such that 

ft = a J rr e jr + aP rs € js and /3 J S = a 3 sr e jr + a J ss e je (10) 

where ar r ,drs etc. are non-zero complex numbers. Moreover, the set Bj — 
{/?r,/?s} is a permutation of B\ — {€j r ,€j s }. As before we write pij = 
(ei\p\ej) = Xij + iyij in the computational basis B\. Then the assignment 
(Bi,Si) fixes the diagonal entries pa and the addition of equations corre- 
sponding to (Bj, Sj) to those of (B±, S\) yields exactly two relations 

Re(a 3 rr a 3 rs p jrjs ) = p> r and Re(ai r ai s p jrjs ) = p> s 

Observe that ((3 r \ (3 S ) = implies that From this it follows 

that the left sides of the two equations are dependent. Since they are con- 
sistent the right sides have similar dependence. That is, the addition of an 
arbitrary assignment (Bj,Sj) G C S' to (Bi,Si) yields exactly one additional 
equation involving Xj r j 3 and Uj r j 3 - Call it the j-incremental equation Sj. 

It follows that the equations for 5f are equivalent to the set of equations 
for (Bi , Si), which fix the diagonal, together with the equations E^, ■ ■ ■ , £f? n +i , 
each of which gives exactly one equation on the real and complex parts of 
some off-diagonal entry. (To get this from the above, note that the equations 
for {(Bi, S^} U X U {(Bj, Sj)} are equivalent to those for {(Bx, Si)} U X U 
{(Bi, Si), (Bj, Sj)}, which are equivalent to {{Bi, Si)} UX U {(Si, Si), Sj}, 
i.e., {(Bi, Si)}U X L){£j}, and apply induction.) Moreover, since rank(^) = 

1 Here the computational basis denotes the standard basis in C n where ei is the vector 
with 0's everywhere except the ith entry which is 1. 
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n 2 , for each off diagonal entry there are exactly two indices j,k in 2 . . . R n 
such that £j and £/- concern that entry. Taking the entry to be that which 
is constrained by £n n+ i, we find that £j and £k and £r„+i concern the same 
entry. But ( S' is equivalent to {B\, S\) with £2, ■ ■ ■ ,£n n , which is consistent, 
and ^ is equivalent to adding to this £r u -\-\, producing an inconsistency. 
The only way this is possible is if £j and £k and £r„+i are inconsistent. 
Since R n > 4 for n > 3, we may find an i? n -subset C S" of §f that contains 
{Bi, Si), (Bj, Sj), (Bk, Sk) and (Bn n+ i, Sr u+ i). But the equations for this set 
then imply the inconsistent set {£j, £k, £n n +i}, so C S" is inconsistent. This 
contradicts the assumption that every i2 n -subset of 5f is consistent. □ 

Observe that the proof of theorem [1] actually provides an algorithm for 
verifying consistency of a given set & of assignments. Thus we start with a 
an arbitrary ^ C & of cardinality R n . If it is inconsistent stop, and declare 
that & is inconsistent. Otherwise find rank(^) = k, say. If k = n 2 then we 
have hermitian matrix p which satisfies the assignments in <S . The problem 
then is to check whether the solution is positive definite. This can be done 
by using the algorithm for Cholesky decomposition [3]. If the algorithm 
succeeds then, of course, p is positive (semi)definite. If it fails then we stop, 
declaring the system inconsistent. Next, suppose k < n 2 . Choose a subset of 
of £f that has maximal rank and let | = s. Add, R n —s new assignments 
and continue the process. If we exhaust & before attaining full rank then 
we have a system of equations with rank r < n 2 . Then the hermitian matrix 
p is determined with n 2 — r free parameters. To determine whether there 
is a positive definite solution in the corresponding parameter space we may 
again follow a Cholesky decomposition type of algorithm but now symbolic. 
In the next section, we present bases in 2 and 3 dimension, inspired by the 
proof of the theorem, to demonstrate that the R n is optimal for dimension 
n. That is, we have determined the consistency number. We conclude this 
section with an application of the Tukey's lemma (equivalent to axiom of 
choice) [5] that yields the following corollary. We use the notations explained 
before. 

Corollary 1. Let & = {(Bj,Sj) \ j £ J} be a collection of quantum prob- 
ability assignments to bases Bj in an n dimensional space, indexed by a set 
J of arbitrary cardinality. Then is consistent if and only if every finite 
subset of cardinality n 2 — n + 1 is consistent. 
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3 Optimality 



In this section we complete the proof that the consistency number R n given 
in Theorem Q] is optimal. We show that there exist quantum probability 
assignments & = {(13%, Sk) \ k = 1, 2, . . . , m > R n } such that every R n — 1- 
subset of & is consistent (hence also every r-subset of & for r < R n is 
consistent) but & is not. For example, in 2 dimensions the bound is 4. 
So we have to show that there exist inconsistent probability assignments 
to 4 bases such that any three of them is consistent. For 3 dimensions the 
bound is 7. We obviously expect the construction of bases and corresponding 
probability assignments to be much harder in 3 and higher dimensions. We 
construct these examples in this section. 

3.1 Dimension=2 

Let the dimension n = 2 and define 1-dimensional projections 



Let the four bases (rather the projections corresponding to the basis 
vectors) be given by 



where / is the 2-dimensional unit matrix. Then the following simultane- 
ous probability assignment is unsatisfiable although any three of them is 
satisfiable. 



In two dimensions we need specify the probability corresponding to only one 
of the outcomes in a projective measurement. It is routine to verify that 
this assignment is satisfiable for any three bases but not for all four. The 
two dimensional case is exceptional as evidenced by the failure of theorems 
of Kochen and Specker and Gleason[6]. 

3.2 Dimension 3 and higher 

In three dimensions the consistency number according to Theorem Q] is 
r'3 < 7. We will construct a projective measurement system consisting of 7 




Bj = {P{,I-P{}, j = l,. ..,4 



p{Pl) = \, p(P?) 



JL p fps\ = 3 and (p4) = 1 
12 » y\ \) g f\ i) 16 
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bases and corresponding probability assignments such that any subset with 
6 bases is consistent but the whole set is not satisfiable showing that the 
consistency number is indeed 7. The bases will be specified by the corre- 
sponding projectors and since each basis is a complete set (the projectors 
add to I, the identity matrix in 3 dimensions) we need specify only two of 
them. Let 

B^iPlPiPl), i = l,...,7 
'0 0^ 

Pi 



p: 



pi 












p 










Vo 





oy 





Po 



Pi 





and P 2 7 



The seven bases given above are particularly simple. We now add another 
base B' = {Qi,Q2\ with more complicated structure. 



Qi 



i 

3 

-7iri/12 

i/3 



\ 3V2 



7iri/12 

1 

2 

g?ri/4 
2V3 



e 7ri/3 \ 

D — 7TZ/4 



2V3 
1 



and Qi 



_6_ 
IT 



,3iri/4 



Now suppose we have the probability assignments 

v(Pj) = p l j and v(Qj) = qj, i = 1,. . . ,7 and j 
If there exists a density matrix 



-3tti/4 

1 

3 

-57TI/12 

^3^ 



1,2. 



e — 7ri/3 

~vT" 

e 5iri/12 

^3^ 
1 



(11) 




cii > and a\ + a2 + 03 = 1 



(12) 
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with measurement probabilities Tr(pPj) = Pj then the following sets of 
equations determine the real parts of Zi. 

a\ = p\ and a 2 = p\ (13a) 

Re{z 3 ) = p\- ( ai + a 2 )/2 = ( fll + a 2 )/2 - p\ (13b) 

Re(z 2 ) = p\ - ( ai + a 3 )/2 = (ax + a 3 )/2 - p| (13c) 

Re(zi) = - (a 2 + o 3 )/2 = (a 2 + a 3 )/2 - ^ (i 3d ) 

Similarly, we have the imaginary parts determined by the bases ,65, Bq, and £7. 



Im(z 3 ) = (ai + a 2 )/2 - p? = p| - (ai + a 2 )/2 (14a) 
Im(z 2 ) = (ai + a 2 )/2 - p\ = p\ - (ax + a 2 )/2 (14b) 
Im(zi) = (ai + o 2 )/2 - p\ = p 2 - (ai + a 2 )/2 (14c) 

Note that the addition of each of the blocks for the bases B 2 , . . . ,Bj to the 
block for basis B\ increases the rank by exactly 1. Finally, for the base 
& = {Qii Q2} we have 

— i7n —iiv i/7T_ 

a\ 0,2 a 3 2Re(z 3 e 12 ) 2Re(z2e 3 ) Re(zie 4 ) 

qi = Y + Y + Y + + 3V2 + ^7T" 

6 fax a 2 2Re(z 3 e tJ r) 2Re(z 2 ef) 2Re(z 1 e^ L )\ 

q2= n{Y + Y + a3+ V6 + V2 + V=s ) 

(15) 

The choice of bases yielding these equations follows a pattern. Each block 
for basis Bi, i > 1 determines exactly one unknown. For example, B3 fixes 
Re(z 3 ) (see ()13b|) ) . Moreover, it is clear that the sets of equations (fT3|) and 
(|14p are independent. So, if two distinct subsets consisting of 6 equations 
each are satisfied by some density matrix then we have already found a 
unique solution. We infer that this set of 7 bases cannot used to show that 
r 3 = 7. Hence we have an additional base B' . Let us then drop one of the 
bases, say B7, and replace it with B' . Let O = {Bx, ■ ■ ■ ,Bq,B'}. The last 
basis B' has the property that if we omit any other basis say Bq from O then 
the remaining system of equations has maximal rank. That is, the resulting 
system of equations (fl~3j) . (|141 ) (omitting (|14cp ) and (fl~5|) is overdetermined. 
We describe below a "procedure" for finding probability assignments to the 
bases in O such that any six of them are consistent but the whole set is not. 
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First, a necessary and sufficient condition for p given in (|12p to be density 
matrix is that the following relations are satisfied. 

3 

det(p) = aia 2 a 3 ai\zi\ 2 + 2Re(ziz 2 z 3 ) > 0, 

aia 2 — \z\\ 2 > 0, 0301 — I ^2 1 2 > 0, a\a2 — \z\\ 2 > and 01,02 > 

These conditions reflect the fact that for a hermitian matrix to be posi- 
tive definite it is necessary and sufficient that all the principal minors have 
non- negative determinant. Now choose a density matrix p such that the in- 
equalities in (fT6|) are strict. We also choose 01,02 > such that a\ + a 2 < 1 
and substitute 03 = 1 — a\ — a 2 . Now treat the Oj, Re(zj) and Im(zi) as 
(real) variables. Then the left side of all the inequalities above are con- 
tinuous functions of these variables. Here we consider p, parametrized by 
{01, 02, Re(zj), Im(zj) : i = 1, 2, 3}, as a member of ]R 8 (R is the field of reals). 
Hence there is an open neighborhood N\ of p such that these inequalities 
hold everywhere in N\. Compute Tr(/oPj) = p 1 -, i = 1, . . . , 6, Tr(pQi) = q\ 
and Tr(pQ 2 ) = q' 2 (note the "prime"). From the set of bases O if we drop 
any base from the set {Bi, . . . , Bq} then the remaining set of equations have 
maximal (real) rank 8. Let Ai denote the invertible matrix of maximal rank 
consisting of coefficients of a set of independent equations corresponding to 
the omission of Bi, i = 2, . . . , 6. 

The invertibility of the matrices Ai, i = 1, . . . , 6 implies that the images 
Ai(Ni) are open sets and hence f] i Ai(Ni) = G is open. The point a' = 
ip\iP\i ■ ■ ■ iP\iP 2 i c Li,Q 2 ) T is i n G as it is the image of p. We choose a point 
q 2 7^ q' 2 such that if the last "coordinate" q' 2 in a' is replaced by q 2 then the 
resulting vector a € G. Since G is open such a choice is always possible. This 
choice of q 2 makes the new assignments for O inconsistent. But restricting 
this assignment to any six bases is still consistent. First, if we omit B' then 
clearly the assignment in rest of the bases is satisfiable by p itself. If we 
omit any other base say Bi then a desired density matrix, say 7, is given by 
A'^ 1 (a). We thus conclude that: 

The probability assignments to the seven bases O given above 
are such that any six of them is consistent with a quantum state 
(density matrix) but the whole set is not. 

Next, it is clear that we can mimic the construction of the bases given 
above for 3 dimensions in any dimension n > 3. Thus, define the bases 
Bo and Bjj, i > j as follows. Let = (0, . . . , 0, 1, . . . , 0) T (ith coordi- 
nated) be the standard basis (the computational basis). Then, Bq = {P® = 
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eiej : i = 1, . . . , n — 1} is the basis consisting of projectors on the standard 
basis. They determine the diagonal elements of the state p. Next, let 

1 < i < j < n and 1 < k < n - 1}, P\ j = -(e* + ej-)(e< + e^) 7 , 
- ej)(ei - ej) T and P l k J = e k e\ k / i,j 

(17) 

These n(n— l)/2 bases determine the real parts of the off-diagonal elements 
of p. These correspond to £>2>£>3 and £>4 in the 3-dim case above. We sim- 
ilarly define the n(n — l)/2 bases for the imaginary parts of off-diagonal 
elements of p. 

1 < j < k < n and 1 < k < n — 1}, P? k = -(ej + iefc)(ej + iek) T , 
iek)(ej - ie k ) T and Pf = e x ej I ^ j, k 

(18) 

As in the 3 dimensional case we replace one of the bases, say B' n _ 1 n , by a 
basis Q. Call the new system of bases Z. The basis Q is chosen so that 
the probability assignments to all the bases in ^ yield an over-determined 
system of equations. This can be achieved by ensuring that the rank the 
system corresponding to — Bij and Z — B'^, VI < i < j < n is maximal 
(= n 2 ). We can then use topological arguments similar to the 3 dimensional 
case to show that every subsystem of Z has a consistent solution but the full 
system Z consisting of R n = n 2 — n + 1 bases does not. Thus, the number 
R n in Theorem Q] is the best possible. We can now state the following. 

Theorem 2. The consistency number of an n- dimensional quantum system 

Combining this theorem with Gleason's theorem [6] we get the following 
theorem. First, recall some definitions needed to state Gleason's famous 
result. Let v be a function on the set V{H) of projections on a finite- 
dimensional Hilbert space H such that 

< v{E) < 1, v(I) = 1 (19a) 
v {E + F) = v(E) + v (F) if EF = (orthogonal projections) (19b) 

Such a function is called frame function. 



P ij - -(e- - 



B' jk = {if : 
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Theorem 3. Let 3? = {Ui \ i G 1} where X is an indexing set and Ui = 
{E\, . . . , E l k .} consists of orthogonal projections: E % -E\ = 5ji on a Hilhert 
space of dimension n > 3. Suppose a real-valued function f on 3? satisfies 
(|19a[) . Further, assume that for every set S C 3f of cardinality n 2 — n + 1 
there exists a frame function v$ such that f(E) = vg(E) for all e G US. 
Then there exists a frame function v such that f = v on L)3f. 

Proof. Gleason's theorem states that any frame function w on a Hilbert 
space of dimension at least 3 is induced by some density operator p: w(E) = 
Tr(pE). This implies that for every subset S of 3f of cardinality n 2 — n + 1 
the function / defines a consistent quantum probability assignment in the 
sense explained in Section [TJ However as the projections Ef may not be 
1-dimensional we cannot apply Theorem [T] directly. In this case we adopt 
the following procedure. Let E be a projection operator of rank k > 1. 
We decompose E = E\ + . . . + E^ where Ei are 1-dimensional projection 
operators with E^Ej = <5jj. We replace each member of Ui by the projectors 
in its 1-dimensional decomposition. Let the resulting set of orthonormal 
projectors be U[. Let 3? = {U-}. It is easy to see that any S' C 3f' 
has consistent probability assignment induced by the original assignment on 
S C 3f . Therefore, from Theorem [T] there is a frame function v (given by 
a density matrix) which yield the same probabilities. The additivity of the 
frame functions now implies that / = v on L)3f. □ 

4 Applications 

In this section we consider some applications of the constructions in the 
previous section. First, we sketch a secret sharing scheme involving 
copies of entangled qubits. A version of the secret sharing problem [7] is 
as follows. A group of k players are to share a secret (represented by a 
number). If any subset of at least r players pool their resources (their shares) 
the secret is revealed otherwise it is not. This is called a (k, r) (threshold) 
secret sharing scheme. Our scheme is probabilistic and we only require that 
any set of r players can discover the secret with high probability whereas 
for less than r players the probability is low. We sketch a scheme using 
entangled qubits. 

Suppose we prepare multiple copies of a composite quantum system of 
dimension N. To each of the k players we send a number of copies of the 
system (or some part thereof) along with instructions for specific measure- 
ments. If we arrange it so that the probability of reconstructing the set 
from any r < k expectation values is high, but negligible for any subset of 
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cardinality r' < r measurements then we have a (A;, r) secret sharing scheme. 
We outline below such a scheme for a system of n qubits. Let N = 2 n , the 
dimension of the system. For simplicity we will only consider a (k, k) secret 
sharing scheme. 

1. A large number of copies of a state randomly chosen from an initial 
set S of Kq states is prepared by the dealer. 

2. The dealer provides the players with m copies each from the original 
ensemble called their shares. 

3. The players are also given instructions about their respective bases 
in which projective measurements are to be performed: the bases are 
from Bij or B'^ given in (fTT|) and (fTKj) . There are N 2 — N such bases. 
Thus each player is actually measuring the real or imaginary part of 
an off-diagonal element of the density matrix. 

4. We assume that the diagonal elements are conveyed to each player 
along with the measurement instruction. These may be used as group 
"password" and/or as a check for interference. 

5. We assume that the set S of states are so chosen that when all but 
one of the (real) parameters characterizing the off-diagonal elements 
of p are fixed there still a large number > K\ of possible states with 
different values of the remaining parameter. We omit the details of 
how this is done in in this sketch. It implies that Kq = 0{K\n 2 ). The 
values of each parameter are separated by a distance > A. 

6. Suppose a is one of the parameters. We see from equations like ([T7]) 
and (I18p that the expectation values yield the probabilities and hence 
the value of a. Using the Chernoff bound (classical) it is seen that if 
ao is the correct value of a and a the calculated value then 

probability (\a — a\ > A) < 0(e~ x /m) 

Therefore if we are aware of the states in S then with high probability 
we can determine the parameter a. 

7. If all the TV 2 — ./V players combine their measurement results then 
(along with the information about diagonal elements) the particular 
state p is determined with high probability. However, with just one 
player missing the probability drops to < \jK\. 
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In this probabilistic protocol we assume that we have secure quantum chan- 
nels with negligible errors. It is possible to devise a more elaborate scheme 
to accommodate insecure channels. Similarly, we can devise a general (k, r) 
secret sharing scheme combining the quantum scheme with a classical one 
involving polynomial evaluation. 

The selection of the set S is a bit more challenging. There are several 
possible approaches however. One is to start with a positive definite matrix 
and then keeping all the values of the parameter vary one of the parameters 
characterizing the off-diagonal elements. For example, let Oy = Re(py) and 
fyj = Im(pij), i < j and suppose we want to vary /3i2 . Using topological 
arguments as before we know that there some neighborhood of 0i2 such that 
for all values of the latter in that neighborhood we get a positive definite 
matrix (fixing the normalization is trivial). We divide the neighborhood 
into segments of appropriate size (fixing A) and pick our values for fii2- 
Alternatively, we could use the fact that for any hermitian matrix A, A 2 is 
positive definite. We vary the parameters of A to achieve our goal. We will 
not go into a detailed analysis here as our primary goal was to demonstrate 
potential applications of the constructions of preceding sections. 

Our protocol has some similarity to the one given in [8] . The difference is 
ours is probabilistic. However, our protocol accommodates a larger number 
of players at the expense of requiring multiple copies of states. Thus using 
4 qubits (N = 16) we can accommodate N 2 — N = 240 players. 

The state p may also be prepared by a purification process [9]. Thus, 
given a density matrix p in an N dimensional Hilbert space we can find 
a pure state in Hm ® Hn, where Hm is an M-dimensional space, such 
that 

P = Tt Hm (\V) <tt|) 

The operation Tr// M is the partial trace with respect to Hm- This method 
eliminates the need for creating a mixed state by random selection. The 
dealer prepares the state l^) and the players get only the components 
(qubits) in H]y. An added advantage of this procedure is that the dealer 
can use the fact that \fy) is entangled with states in Hm to ensure that the 
players are using their assigned measurements. 

The bases given in the equations (|17p and (|18p could be used for quantum 
state tomography. In this case, we do not require maximal outcomes in all 
the bases. Thus, assuming the diagonal elements of the density matrix have 
been estimated we need only measurements with three alternatives. For 
example, if we want to estimate Re(/9i2) then using the notation in (I17p the 
orthogonal projections P 1 , P 2 12 and M 12 = J2k^i 2 -^fc 2 f° rm a complete set 
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and provide the three outcomes. Note also that projections P- J and Pj 3 can 
be obtained by repeated applications of Hadamard and control gates. 

5 Discussion 

In this work we stated and solved a problem of consistent probability as- 
signments for maximal projective measurements, that is, the number of pro- 
jection operators in each measurement equals the dimension. This can be 
relaxed. For example, in the measurement bases in Section 13.21 it is imma- 
terial for the basis Bi, i > whether we take the complete basis or three 
projectors — two orthogonal projectors on a "plane" spanned by two vec- 
tors from Bq and one projector orthogonal to the plane. The problem of 
consistency number r n arises in the firs place because we have no a priori 
knowledge about the independence of the projection operators correspond- 
ing to the bases. We can visualize such a situation when different observes 
have no initial communication and can perform only local measurements. 
Further, it is not easy to define general procedures for constructing such 
independent bases. One notable exception is the explicit recipe for mutually 
unbiased bases (MUB). MUBs were first introduced by Schwinger [10] in 
low dimension and later extended to prime and prime power dimensions by 
Ivanovic [11] and Wootters and Fields [12] respectively. The problem is this 
recipe does not work when the dimension has two distinct prime factors. 

The next logical question would be to consider the problem of consistent 
probability assignment for more general measurement schemes, in particu- 
lar, for local measurements of entangled states. The most general problem 
regarding consistency number would be the following. What is the consis- 
tency number of a collection of measurements (positive operator valued in 
general) given some prior information about the state? In this format we 
have to satisfy some extra constraints. For example, we may have informa- 
tion that the unknown state is pure. The constraints would be nonlinear in 
general. More interesting protocols would result from these investigations. 
We hope to investigate these questions in future. 

References 

[1] van der Meyden R and Patra M 2003 In Proc. Computer Science Logic 
and 8th Kurt G'bdel Colloquium, Vienna, 2003. Springer- Verlag 

[2] Bell J S 1966 Rev. Mod. Phys. 38 447 



18 



[3] Kochen S and Specker E P 1967 Jour. Math, and Mechanics 17 59 

[4] Stewart G W 1998 Matrix Algorithms Vol. 1, SIAM, Philadelphia p. 188 

[5] Kelley J L 1975 General Topology, Springer p. 33 

[6] Gleason A M 1957 Jour. Math, and Mechanics 6 885 

[7] Shamir A 1979 Comm. ACM 22 612 

[8] Hillery M, Buzek V and Berthiaume A 1999 Phys. Rev. A 59 1829 

[9] Nielsen M A and Chuang I L 2000 Quantum computation and infor- 
mation, Cambridge 

[10] Schwinger J 1960 Proc. Natl. Acad. Sci. USA 46 570-579 

[11] Ivanovic I D 1981 J. Phys. A 14, 3241 

[12] Wootters W K and Fields B D 1989 Annals of Physics 191, 363 



19 



